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ABSTRACT 

We examine the metallicity distribution function (MDF) and fraction of carbon-enhanced metal-poor 
(CEMP) stars in a sample that includes 86 stars with [Fe/H] < -3.0, based on high-resolution, high-S/N spec- 
troscopy, of which some 32 objects lie below [Fe/H] = -3.5. After accounting for the completeness function, 
the "corrected" MDF does not exhibit the sudden drop at [Fe/H] = -3.6 that was found in recent samples of 
dwarfs and giants from the Hamburg/ESO survey. Rather, the MDF decreases smoothly down to [Fe/H] = -4. 1 . 
Similar results are obtained from the "raw" MDF. We find the f ract ion of CEM P objects below [Fe/H] = -3.0 is 
23 ± 6% and 32 ± 8% when adopting the lBeers & Christliebl and lAoki et ail CEMP defin itions, respectively. 
The former value is in fair agreement with some previous measurements, which adopt the iBeers & Christliebl 
criterion. 

Subject headings: Cosmology: Early Universe, Galaxy: Formation, Galaxy: Halo, Nucleosynthesis, Abun- 
dances, Stars: Abundances 



1. INTRODUCTION 

Metal-poor stars provide critical information on the ear- 
liest phases of Galactic form ation (see e.g., the rev iews by 
IBeers & Christliebl 120051 and IFrebel & Norrisl l20Tlb . Their 
chemical abundances shed light upon the nature of the first 
stars to have formed in the Universe, and the nucleosynthesis 
which seeded all subsequent generations of stars. 

This is the third paper in our series, which focuses upon 
the discovery of, and high-resolution, high signal-to-noise ra- 
tio (S/N) spectroscopic analysis of, the most metal-poor stars. 
Here we explore two key issues: the metallicity distribution 
function (MDF) and the fraction of carbon-enhanced metal- 
poor (CEMP0 stars at lowest metallicities. 

1 This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 

2 Some of the data presented herein were obtained at the W. M. Keck 
Observatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and the 
National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W. M. Keck 
Foundation. 

3 Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile (proposal 
281.D-5015). 

4 Research School of Astronomy and Astrophysics, The Australian Na- 
tional University, Weston, ACT 2611, Australia; yong@mso.anu.edu.au, 
jen@mso.anu.edu.au, bessell@mso.anu.edu.au, martin@mso.anu.edu.au 

5 Zentrum fur Astronomie der Universitat Heidelberg, Lan- 
dessternwarte, Konigstuhl 12, D-69117 Heidelberg, Germany; 
n.christlieb @ lsw.uni-heidelberg.de 

6 Max-Planck Institute for Astrophysics, Karl-Schwarzschild Str. 1 , 
85741, Garching, Germany 

7 National Optical Astronomy Observatory, Tucson, AZ 85719 

8 Department of Physics & Astronomy and JINA: Joint Institute for 
Nuclear Astrophysics, Michigan State University, E. Lansing, MI 48824, 
USA; beers@pa.msu.edu 

9 Department of Physics and Astronomy, Uppsala University, Box 515, 
75120 Uppsala, Sweden; paul.barklem@physics.uu.se 

10 Massachusetts Institute of Technology, Kavli Institute for As- 
trophysics and Space Research, Cambridge, MA 02139, USA; 
afrebel @ mit.edu 

11 Centre for Astrophysics Research, School of Physics, Astronomy 
& Mathematics, University of Hertfordshire, College Lane, Hatfield, 
Hertfordshire, ALIO 9AB, UK; s.g.ryan@herts.ac.uk 

12 Initially defined as stars with [C/Fe] > +1.0 and [Fe/H] < -2.0 



Any model purporting to explain the formation and evolu- 
tion of our Galaxy must be able to reproduce the observed 
MDF. The ingredients of such models include the initial mass 
function (IMF), nucleosynthetic yields, and inflow or outflow 
of gas. Observations of the MDF can constrain these initial 
condition s and p hysical processes. Since the early work by 
iHartwickl dl976l) . measurements of the MDF involve increas- 
ing numbers of stars with more accurate metallicity me asure- 
ments (see e.g. jLaird etai]|1988l:lR"van & Norrisll 19911) . One 
of the basic predictions of Hartwick's Simple Model of Galac- 
tic Chemical Enrichment is that the number of stars having 
abundance less than a given metallicity should decrease by 
a factor of ten for each factor of ten decrease in metallicit^Q 
iNorrisI (1 1 9991) presented observational support for this sugges- 
tion, down to [Fe/H] ~ -4.0 ? below which i t appe ared to be no 
longer valid. More recently. lSchorck et ail (120091) and lLi et ail 
(120101) presented MDFs of the Galactic halo using 1638 gi- 
ant and 617 d warf stars, respectively , from the Hamburg/ESO 
Survey (HES: IWisotzki et aflfl996l) . Below [Fe/H] = -2.5, 
the MDFs for dwarfs and giants were in excellent agreement. 
A prominent feature of both MDFs was the apparent lack of 
stars more metal-poor than [Fe/H] = -3.6. While a handful of 
such stars are known, the sharp cutoff in the MDF has impor- 
tant implications for the critical metall icity above which low - 
mass star formation is possible (e.g., ISalvadori et al.ll2007l) . 
More detailed studies of the MDF, and in particular the low- 
metallicity tail, are necessary to confirm and constrain the 
star fo rmation modes of the first stars (e.g.. lBromm & Larsonl 
2004). 

The HK survey (IBeers et al.l 19851 1 19921) revealed that there 
is a large fraction of metal-poor stars with unusually strong 
CH G-bands indicating high C abundances. With the addition 
of numerous metal-poor stars found in the HES, the CEMP 
fraction at low metallicity has been c onfirmed and qua nti- 
fied, wit h estimates ranging f rom 9% (IFrebel et al .1120061) to 
> 21% dLucatello et al.1 12006). These numbers are consid- 

IBeers & Christlierj|2005l) . 

13 While a number of c hemical evolution mod el s (e.g., iKobavashi et al.1 
20061 IKarlssonl 12 006. Sal vadori et aT] [2007T. IPrantzosI \205S. and 
Cescutti & Chiappini 20101) have improved upon the one-zone closed- 
box Hartwick model, the general behavior remains largely unchanged. 
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erably larger than the fraction of C-rich objects at higher 
metallicity, the so-called CH and Ba stars, which account for 
only ~ 1 % of the population. The fraction is even larger at 
lowest metallicity: below [Fe/H] < -4. 5, 75% of the four 
known stars belon g to the CEMP class dNorris et al.l 120071 
ICaffau et al.ll2.Ql lb . To explain these large fractions, several 
studies argue that adjustments to the IMF are necessary (e.g ., 
Lucatello et al]|200l iKomiva et al.ll2007t llzzard et aljl2009l 
Carollo et al.l (120121) offer an alternative interpretation for the 
increase of the CEMP fraction they observe in the range -3.0 
< [Fe/H] < -1.5 in terms of a dependence of CEMP fraction 
on height above the Galactic plane. In their most metal-poor 
bin at [Fe/H] ~ -2.7, they report C-rich fractions of 20% and 
30% for their inner- and outer-halo components, respectively 
(see their Figure 15). In their view, this can be accounted for 
by the presence of different carbon-production mechanisms 
(some not involving the presence of AGB nucleosynthesis) 
that have operated in the inner- and outer-halo populations. 

An understanding of the CEMP stars is complicated by 
the fact that they do n ot form a homogeneous group: 
iBeers & Christliebl (12005b define several CEMP subclasses 
(all of which have [C/Fe] > +1.0) as follows: (i) CEMP-r - 
[Eu/Fe] > +1.0; (ii) CEMP-s - [Ba/Fe] > +1.0 and [Ba/Eu] > 
+ 0.5; (hi) CEM P r/s - 0.0 < [Ba/Eu] < +0.5; and CEMP-no 
- [Ba/Fe] < 0.0. 1 Aokil (120101) shows that below [Fe/H] = -3.0, 
the CEMP stars are principally (90%) CEMP-no stars, while 
for [Fe/H] > -3.0, the CEMP-s class predominates. These 
differences lie outside the scope of the present paper. Here 
we seek to constrain only the fraction of CEMP stars at lowest 
abundance, [Fe/H] < -3.0, and to compare the results with the 
fractions determined at higher abundances. In Paper IV (Nor- 
ris et al. 2012b) we shall address the nature of the CEMP-no 
stars, which comprise the large majority of CEMP stars in our 
extremely metal-poor sample. 

2. OBSERVATIONS AND ANALYSIS 

In Norris et al. (2012a; Paper I), we presented high- 
resolution spectroscopic observations of 38 extremely metal- 
poor stars ([Fe/H] < -3.0; 34 newly discovered), obtained us- 
ing the Keck, Magellan, and VLT telescopes, including the 
discovery and sample selection, equivalent-width measure- 
ments, radial velocities, and line list. In Paper I, we also de- 
scribed the temperature scale, which consists of spectropho- 
tometry and Balmer-line analysis. In addition to the 38 pro- 
gram stars, we se lected 207 stars from the SAGA database 
(ISuda et al.l 120081) (queried on 2 Feb 2010), and performed 
a homogeneous re-analysis of this literature sample. All 
stars were analyzed using the NEWODF gr id of ATLAS9 
model atmospheres (Ca stelli & Kuruczll2003l) . and t he 2011 
versio n of the stellar line-analysis program MOOG dSnedenl 
119731) . which includes pro per treatment of continuum scat- 
tering (ISobeck et al.l 120171) . They thus have effective tem- 
peratures, surface gravitie s, microturbulent ve locities, log gf 
values, solar abundances (lAsplund et al.ll2009l) . and therefore 
metallicities, [Fe/H], all on the same scale. 

The literature sample was reduced from 207 to 152 stars 
by (i) discarding stars with fewer than 14 Fe I lines (the mini- 
mum number of Fe I lines measured in our program stars), (ii) 
removing literature stars included in the program-star sample, 
and (iii) averaging the results of stars having multiple analy- 
ses into a single set of abundances. Thus, the final combined 
sample consists of 190 stars (38 program stars and 152 litera- 
ture stars). Full details regarding the analysis are presented in 
Yong et al. (2012, Paper II). 



3. RESULTS 
3.1. Selection Biases 

In Table Q] we present data, based on our high-resolution 
analyses, for the 8<£3 stars m our collective sample that have 
[Fe/H] < -3.0; of these, 32 have [Fe/H] < -3.5, while there 
are nine with [Fe/H] < -4.0. We stress again that these metal- 
licities are on our homogeneous system of T e s, logg, logg/ 
values, and solar abundances. These are the most metal-poor 
stars known in our Galaxy, and allow us to address below the 
key issues of the MDF and CEMP fraction. 

Before continuing, we comment on the completeness func- 
tion and selection biases of the sample. The HES is complete 
for metalliciti es below [Fe/H] = -3.0 (ISchorck et alJ 120091: 
iLi et alj|2010|). To e stimat e the completeness. ISchorck et all 
(120091) and lLietal.1 (120101) used the Simple Model to gener- 
ate a metallicity distribution function and then applied their 
selection criteria to obtain the MDF wh ich would have been 
observed in the HES (see Secti on 6 in ISchorck et al.l I2009L 
and Section 3.4 in lLi et a02010l for further details). From Pa- 
per I, we can compute the completeness function for the ~ 30 
HES candidates having high-resolution, high-S/N spectra dis- 
covered in that work. First, we use a linear transformation 
to place the medium-resolution metallicities, [Fe/H]K, from 
Paper I onto the high-resolution abundance scale, [Fe/H]. We 
then compare the number of HES stars observed at high reso- 
lution with the total number of HES stars observed at medium 
resolution, and from which the stars observed at high resolu- 
tion were selected, as a function of [Fe/H]. We use this ratio 
to correct the MDFs in the following subsection. In a simi- 
lar manner, we are able to determine the completeness func- 
tion for the ~ 50 HK-survey stars in our extended sample, by 
using material in the medium-resolution HK database main- 
tained by T. C. B. 

3.2. The Metallicity Distribution Function (MDF) 

Our MDFs are presented in Figure [TIPl where in the left 
panels the scale of the ordinate is linear and for those on the 
right it is logarithmic. The two uppermost panels each contain 
MDFs constructed from the raw data for the 38 program stars 
and the total sample of 190 objects. We use generalized his- 
tograms, in which each data point is replaced by a Gaussian 
of width a = O.3C0 dex. The Gaussians are then summed to 
produce a realistically smoothed histogram. 

Construction of our smoothed MDF includes uncertainties, 
which we estimate in the following manner using Monte Carlo 
simulations. We replaced each data point, [Fe/H], with a ran- 
dom number drawn from a normal distribution of width 0.15 
dex, centered at the [Fe/H] of the given data point. We re- 
peated this process for each data point in our collective sam- 
ple of 190 stars, and a generalized histogram was constructed 

14 For nine program stars, we could not determine whether they were 
dwarfs or subgiants. For the subset of those stars included in this paper, 
we present the results for both cases in Table[T] In all figures, unless noted 
otherwise, we adopt the average [Fe/H] and [X/Fe] from the dwarf and sub- 
giant analyses for these stars. For the nine objects, the average differences are 
([Fe/H] dwarf - [Fe/H] subgiant ) = 0.02 ± 0.01 dex (a = 0.03) and ([X/Fe] dwarf - 
[X/Fe] subgiant ) = 0.05 ± 0.02 dex (a = 0. 17), where X refers to the 14 species 
(from Na to Ba) measured in Paper II. For C, while the differences are larger, 
[C/Fe] sl , b giant - [C/Fe] dwarf = 0.23 ± 0.05 dex (cr = 0.13), the CEMP classifi- 
cations do not depend on whether we adopt the dwarf or subgiant value. 

15 All figures were generated using the full sample, presented in Paper II. 

16 We regard our typical uncertainty in [Fe/H] to be 0.15 dex, rather than 
0.30 dex. Given our still relatively limited sample size, using a = 0.15 dex 
produces spurious structure in our MDF. None of our conclusions depend 
upon our choice of a in constructing the MDF. 
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TABLE 1 

Stellar Parameters and Carbon Abundance 



Star 


RA2()00 a DEC2000 a 


'cIT 


l°gg 


6 


[M/H] m0 del 


[Fe/H] delived 


[C/Fe] b 


C-rich c 


Source 








(K) 


(cgs) 


(kms -1 ) 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


CS29527-015 


00 29 10.7 - 


19 10 07.2 


6577 


3.89 


1.9 


-3.3 


-3.32 


1.18 


1 
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CS30339-069 


00 30 16.0 - 


35 56 51.2 


6326 


3.79 


1.4 


-3.1 


-3.05 


0.56 
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CS29497-034 


00 41 39.8 - 


26 18 54.4 


4983 


1.96 


2.0 


-3.0 


-3.00 


2.72 


1 


4 


HD4306 


00 45 27.2 -09 32 39.9 


4854 


1.61 


1.6 


-3.1 


-3.04 


0.11 





12 


CD-38 245 


00 46 36.2 - 


37 39 33.5 


4857 


1.54 


2.2 


-4.2 


-4.15 


< -0.33 
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HE0049-3948 


00 52 13.4 - 


39 32 36.9 


6466 


3.78 


0.8 


-3.7 


-3.68 


<1.81 





1 


HE0057-5959 


00 59 54.0 - 


59 43 29.9 


5257 


2.65 


1.5 


-4.1 


-4.08 


0.86 


1 


1 


HE0102-1213 


01 05 28.0 - 


11 5731.1 


6100 


3.65 


1.5 


-3.3 


-3.28 


<1.31 





1 


CS22183-031 


01 09 05.1 -04 43 21.1 


5202 


2.54 


1.1 


-3.2 


-3.17 


0.42 





12 


HE0107-5240 


01 09 29.2 - 


52 24 34.2 


5100 


2.20 


2.2 


-5.3 


-5.54 


3.85 


1 
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References. 
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Aoki et alj j 
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et alj (20071 120091); 6 


= Carretta et al. (2002); Cohen et al. 
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2004); Spite etui. 2005); 8 = Christlieb et al. i 2004 
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. 11 =IFrebel et alj 


2007b); 12 = Honda et al. (2004); 13 = Lai et al. i 2008 ); 14 = 


INorriset all 12001 


); 15=INorrisetalJI200" 





Note. Table 1 is published in its entirety in the electronic edition of The Astrophysical Journal, 
its form and content. 

a Coordinates are on the 2MASS system l Skrutskie et al. 2006) 

For literature stars, [C/Fe] is the (average) value from the reference(s). 
c We adopt the lAoki et aT] 120071) CEMP definition. 

This analysis assumes the star is a dwarf. 
e This analysis assumes the star is a subgiant. 
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for this new sample. We repeated this process for 10,000 new 
random samples, producing a generalized histogram for each 
new random sample. At a given [Fe/H], we then have a dis- 
tribution of some 10,000 values, one for each MDF. We mea- 
sured the FWHM of this distribution, and adopt this value as 
an estimate of the uncertainty in our MDF at a given [Fe/H]. In 
Figure 1(c), we plot the fractional uncertainty, where a value 
of 0.2 represents a 20% uncertainty in the value of the MDF. 
The relative uncertainty reaches 50% near [Fe/H] = -4.2, and 
becomes rapidly larger at lower metallicities, indicating that 
the sample size loses much statistical significance below this 
value. 

We also constructed a regular histogram to com- 
pare with the smoo th ed M DF. We employed the 
IShimazaki & Shinomotol J2007I) algorithm to determine 
the optimal bin width (0.272 dex) for the full 190 star sample. 
As expected, both histograms exhibit a similar behavior. 

We corrected the "program star MDF" using the HES 
completeness function described above in Section 13.11 (here 
shown together with the HK completeness function in Figure 
1(d)), leaving the "literature sample MDF" unchanged. These 
MDFs are presented in Figure Q] (panels e-f). We also cor- 
rected the full MDF (i.e., "program star + literature sample" 
MDF) using the HES completeness function, and plot both 
corrected MDFs in Figure \T\ (panels g-h). While the selec- 
tion biases associated with the discovery of the stars in the 
SAGA database are not explicit, almost half of the 86 stars 
(42) in Table Q] carry HK-survey names, while most others 
(36) have HES-survey nomenclature. It is clear that the major- 
ity of stars in Table 1 have been found in those low-resolution 
spectroscopic surveys, and thus inherit the spectroscopic- and 
volume-selection biases of those works, plus additional biases 
imposed in later follow-up with medium- and high-resolution 
spectroscopy. Many of the HK-survey stars would also have 
been recovered in the HES survey, but were not renamed. 
Consequently, using the HES completeness function should 
be a reasonable step. Given the clear similarity between 
the HES and HK completeness functions below [Fe/H] = 
-3.3, the corrected MDF would be essentially identical in this 



metallicity regime had we used the HK completeness func- 
tion. 

We use a two-sample Kolmogorov-Smirnov (KS) test to 
compare the MDFs for dwarfs (logg > 3.5) and giants (logg 
< 3.5). The null hypothesis is that the dwarf and giant MDFs 
are drawn from the same distribution. For [Fe/H] < -3.0, the 
two-sample KS test yields a probability of 0.601 (D = 0.167) 
that the dwarf and giant MDFs are drawn from the same dis- 
tributiorQ A similar test for [Fe/H] < -3.5 yields a prob- 
ability of 0.915 (D = 0.200) that the dwarf and giant MDFs 
are drawn from the same distribution. Therefore, the null hy- 
pothesis that the giants and dwarfs are drawn from the same 
population cannot be rejected at the .10 level of sig nificance, 
the least stringent level in Table M of lSiegell d 19561) . 

In Figure ELTa), we overplot the raw MDF from 
ISchorck et al.1 (120091) (usi ng the values in their T able 3 ). Com- 
paring our sample with ISchorck et al.l (120091) and iLi et alJ 
(120101) (made available by N. C), we find 12 stars in com- 
mon. For these 12 stars, there are some 18 [Fe/H] measure- 
ments that can be compared. For the nine program stars for 
which we conducted dwarf and subgiant analyses, we treat 
both [Fe/H] values as independent measurements for the pur- 
poses of this compa rison . Our metallic ities differ from the 
ISchorck et aD (120091) and lLi etalJ (120101) values by -0.26 ± 
0. 06 dex (cr = 2 7 dex ), and so we shift the raw HES MDF 
of ISchorck et alj ( 120091) by -0.26 dex in Figure Q] and scale 
it to match our MDF at [Fe/H] = -3.5. Below [Fe/H] = -3.5, 
we find a large fraction of stars relative to the HES MDF. In 
FigureQIf), both the program star and literature sample MDFs 
have a slope close to 1.0, consistent with the Hartwick Sim- 
ple Model, down to the shoulder at [Fe/H] ~ —4.1, when the 
finite sample begins to run out of stars (which are necessar- 
ily counted in integers). This corresponds to the metallicity at 
which the fractional error (FigureHlc)) increases rapidly and, 
as noted above, the finite sample size loses much statistical 
significance. 

Therefore, taken at face value, and bearing in mind the bi- 

17 The dwarf and giant MDFs for [Fe/H] < -3.0 may be seen in Figure|6] 
panels (d) and (e), respectively, which we shall discuss in what follows. 
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1 5 1 (a) Raw MDFs 

All stors 
(N = 190) 
Program stors 
(N=38) 
0.5 h HES 




-5 -4 
[Fe/H] 

FIG. 1. — Generalized histograms showing the MDF (linear (left) a nd logarithmic (right) scales). The full sample (solid black line) and program stars (red 
histogram) are shown. The green dashed line is the raw HES MDF from Schorck et al. 1 2009), shifted by -0.29 dex and scaled to match our MDF at [Fe/H] = — 3.5. 
Panel (a) includes a regular histogram (dotted line) employing the Shimazaki & Shinomoto 1 2007) optimal bin width algorithm. Panel (c) shows the fractional 
uncertainty in our MDF (e.g, a fractional uncertainty of 0.2 represents an error of 20% of the MDF value.) Panel (d) shows the HES and HK completeness 
functions. The HES completeness function is applied to the MDFs shown in panels (e,f,g,h). 



ases, the apparent cutoff in the HES MDF at [Fe/H] = -3.6 
is not confirmed by our data. We identify 13 HES stars in 
our sample that have [Fe/H] < -3.7 (of w hic h four are con - 
tained in the work of ISchorck et al.H2009l and lLi et al.ll2010t) . 
We speculate that (i) stars in our sample having [Fe/H] < - 
3.7 were rejected as having strong G-bands (GF0 > 6 A), 
an d/or (ii) our abundanc e sc ale differs from that adopted in 
the ISchorck et al.1 (120091) and lLi et all (120101) analyses. 

Regarding point (i), none of our objects has GP > 6A. In 
particular, we note that the three most Fe-poor HES stars, all 
of which have large [C/Fe] ratios, are not rejected by this cri- 
terion. Concerning point (ii), Figure [2] shows the metallic- 



18 This is the Beers et a I] (1993) index that measures the strength of the 
4300A CH molecular features. 



ity difference A [Fe/H] = [Fe/H] (high resolution: this study) 
— [Fe /Hlic (medium resolution: ISchorck et alj|2009t iLi et al.l 
|20TO|) versus r eff , logg, [Fe/H], E(fl-V), and GP. In each 
panel, we plot the linear least squares fit to the data, and show 
the formal slope and uncertainty as well as the dispersion 
about the slope. In all cases, the dispersion about the slope is 
compatible with the value expected based on the convolution 
of the errors, a (c ombined) = 0.25 de x assuming a (this study) 
= 0.15 dex and a (iBeers et al.ll 19991) = 0.20 dex. The correla- 
tion between A [Fe/H] versus [Fe/H] is significant at the 2-cr 
level (although we caution that the errors on these two quan- 
tities are correlated), indicating that as one moves to lower 
metallicity, the [Fe/H] values from our high-resolution anal- 
ysis are lower than those based on medium-resolution spec- 
tra. Such a correlation would help, in part, to explain why 
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we do not find a cutoff in the MDF. Possible explanations for 
this correlation include systematic differences in the analyses, 
interstellar Ca absorption lines, and/or CH molecular stellar 
lines in the region of the Ca II K line. Further insight into 
the differences from high-resolution and medium-resolution 
spectra await larger comparison samples. (We also note that 
Figure |2] includes 3-a correlations between A[Fe/H] and T e ff 
(panel a) and A[Fe/H] and logg (panel b).) For completeness, 
we note that linear regression analysis shows that the best fit 
to A[Fe/H] (high resolution - medium resolution) is -1.825 
+ 2.442 x 10" 4 x 7k + 0. 15 1 x logg +0. 160 x [Fe/H] hlg h resolution 
- 0.401 xE(B-V) + 0.337xGP. The dispersion about this 
fit is 0.17 dex, and the uncertainties in the coefficients are 
2.986 xlO" 4 , 0.147, 0.189, 6.637, and 0.295 for T e{{ , logg, 
[Fe/H], E(B-V), and GP, respectively. 

In Figure [3] we compare the raw and corrected MDFs 
with several model predictions, scaled to match our MDFs at 
[Fe/H] = -3.5. The rationale for choosing this normalization 
is that (i) in this metallicity regime we expect that our sam- 
ple includes the vast majority of stars currently known, albeit 
with selection biases, and (ii) we hope to provide a more de- 
tailed consideration of the MDF at the lowest observed [Fe/H] 
values. 

All predictions, except the Kobavas hi et al.l ([2006) "out- 
flow" model , provide a reasonable fi t to the raw and corrected 
MDFs. The iKobavashi et al.l d2006l) "infall" model provides 
a superior fit to our MDF than their "outflow" model (which 
overpredicts the number of metal-poor stars). The "outflow" 
model contains (i) outflow, (ii) no infall, and (iii) a low star- 
formation efficiency, while the "infall" model contains (i) no 
outflo w, (ii) infall, and (iii) a much lower star-formation effi- 
ciency. iPrantzosI d2008l) adopts a hierarchical merging frame- 
work in which the halo is formed from sub-halos, with a dis- 
tribution in stellar mass, and with the MDF of each sub-halo 
based on Local Group dwarf satellite galaxies. Both Prantzos 
models ("outflow" onl y and "outflow+ i nfall" ) provide equally 
good fits to our MDF. ISalvadori et al.l d2007) provide predic- 
tions for different critical metallicities, Z cr , and their Z cl = 
10~ 4 Z© and Z cl = models both provide reasonable fits to our 
MDF. The raw and corrected MDFs indicate that the critical 
metallicity, above which low-mass star formati on is possible, 
is well bel ow Z cr = 10~ 3 4 Z ^, in contrast to the lSchorck et al.l 
(120091) and lLi et all (120101) MDFs. In addition to the spectro- 
scopic selection biases noted earlier, we need to be mindful 
of possible volume-selection biases, and that the real Galactic 
MDF at low metallicities could be significantly different from 
the one presented in this paper. Still larger, deeper samples, 
the biases and completeness of which are better understood, 
are necessary to obtain this MDF. 

3.2.1. On the nature of the MDF 

We now explore four aspects of our MDF analysis: 
(1) choice of a lower-metallicity cutoff versus a higher- 
metallicity cutoff, (2) usage of a regular histogram versus a 
generalized histogram, (3) adoption of a linear versus a loga- 
rithmic scale, and (4) inclusion of elements in addition to Fe 
when defining the metallicity. 

Lower-metallicity cutoff versus higher-metallicity cutoff 

In order to explore the first aspect, we adopt the (one- 
zone , closed-box) S imple Model of Galactic chemical evolu- 
tion dSchmidlll 19631: ISearle & Sargentll 19721: iPagel & Patchettl 



1 1 9751: lHartwicklfr976l) . and create two MDFs, from which we 
remove all stars below [Fe/H] = -4.5 (lower-metallicity cut- 
off) and -4.0 (higher-metallicity cutoff, sometimes referred to 
as a "sharp cutoff"). Both are populated with stars on a regular 
grid of step size 0.05 dex, and normalized such that they have 
1000 stars below [Fe/H] =-3.0,i.e., some 12 times larger than 
our 86 star sample in Table Q] For the lower-metallicity cut- 
off (MDF1), there are four stars in the lowest metallicity bin, 
[Fe/H] = -4.5, while for the higher-metallicity cutoff (MDF2), 
there are 20 stars in the lowest metallicity bin, [Fe/H] = -4.0. 
The two MDFs are shown in Figure|4] In the upper panel, one 
sees that when overplotted on the full metallicity range, -5.0 
< [Fe/H] < 0.0, they are indistinguishable. When considering 
only the regime below [Fe/H] = -3.0 (Figure 4 panels (b) and 
(c)), however, the difference in the two MDFs is clear. 

Regular histogram versus generalized histogram 

Panels (b) and (c) of Figure [4] show regular histograms for 
the two MDFs, while panels (f) and (g) show generalized his- 
tograms. As expected, the generalized histogram smooths out 
the data along the abscissa. Given the numbers of stars in the 
lowest metallicity bins, the lower-metallicity cutoff MDF may 
appear to have an "extended tail," when represented in gener- 
alized histogram format, but in reality, both MDFs now have 
an additional tail. 

Linear versus logarithmic scale 

Panels (b,c,f,g) and (d,e,h,i) of Figure|4]have linear and log- 
arithmic scales respectively. Panels (d) and (e) (regular his- 
tograms) and panels (h) and (i) (generalized histograms) ex- 
hibit rather similar trends. When using a logarithmic scale, 
it is easier to discern where the MDF cuts off, (as every fi- 
nite sample, observed or simulated, must). The generalized 
histogram replaces each datum with a Gaussian function, and 
taking the logarithm of this yields an inverted quadratic func- 
tion; i.e., each datum contributes an inverted quadratic func- 
tion to the log panel. In Figure SJh), the last Monte Carlo 
datum at [Fe/H] = -4.5 gives rise to the quadratic roll-off at 
[Fe/H] < -4.5, and in Figure Hi) the last Monte Carlo da- 
tum at [Fe/H] = -4.0 gives rise to the roll-off at [Fe/H] < 
-4.0. This roll-off meets the populated part of the MDF at 
a "shoulder", above which the MDF rises with a slope of 
1.0, due to the adoption of the Simple Model. The location 
of the shoulder indicates the metallicity at which either the 
finite sample size becomes too small to populate the MDF, 
as in this simulation, or the MDF genuinely departs from 
the Galactic chemical evolution model pertaining at higher 
m etallicity, as would be the case in the scenarios envisaged 
by ISalvadori et al.l d2007l) and others discussed in connection 
with Figure|3] The fact that the shoulder in our observed MDF 
(e.g., Figures Hlf), 1(h) or Figure |3J, determined from high- 
resolution spectroscopic analyses, is located at [Fe/H] = -4.1 
or -4.2, and attains a slope close to 1.0 at higher metallic- 
ity, gives us the confidence that the MDF does not exhibit 
a sharp drop at [Fe/H] = -3.6, nor indeed in the metallicity 
range down to [Fe/H] = -4.1. 

Inclusion of elements in addition to Fe in the "metallicity" 

Strictly defined, metallicity (Z) includes all elements heav- 
ier than helium, although in practice Fe is widely adopted as 
the canonical measure of stellar metallicity. Therefore, the 
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FIG. 2. — The difference in metallicity (This Study - Literature) between our analysis and those of Schorck et al. (2009) (crosses) and Li et al. 1 2010) (circles) 
vs. (a) r e ff, (b) logg, (c) [Fe/H], (d) E(B-V), and (e) GP, where the abscissa values in panels (a), (b), and (c) were obtained from the high-resolution analysis. 
In each panel, we plot the lineal' least squares fit to the data, and show the slope, uncertainty, and dispersion about the slope. In this figure, we include both the 
dwarf and subgiant [Fe/H] measurements for those program stars with multiple analyses (see Section[2]for details). 

MDF discussed thus far is really the Fe distribution func- 
tion. For the Sun, the seven most abundant met als, in decreas- 
ing o rder, are O, C, Ne, N, Mg, Si, and Fe (lAsplund et al.l 
120091) . Therefore, in order to explore this fourth aspect of 
our discussion, the behavior of the MDF when including ad- 
ditional elements, we arbitrarily define Z to consist of C, N, 
Mg, Si, and Fe. (Of the 86 stars with [Fe/H] < -3.0, there 
are measurements of C, N, Mg, and Si for 54, 36, 81, and 
36 stars respectively.) We compute Z for each star, only con- 
sidering the set of elements with measurements; that is, we 



ignore those elements not measured in a given star. In Fig- 
ure^ we plot [Z/H] versus [Fe/H], including all stars in our 
sample (N = 190). Panels (b) and (c) in Figure [5] show the 
MDFs for [Fe/H] and [Z/H], respectively, in regular and gen- 
era lized histogram format. The reg ular histogram again uses 
the IShimazaki & Shinomotol d2007l) optimal bin width algo- 
rithm (0.278 dex for [Z/H]). We note that the two MDFs ex- 
hibit a similar behavior. Indeed, the [Fe/H] and [Z/H] MDFs 
have almost identical average gradients over the plotted range. 
The construction of the [Z/H] MDF, based on large samples of 
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stars having O and C measurements, would be of great inter- 
est given the postulated importance of thes e elements for low- 
mass star formation in the early Universe (iBromm & Larsonl 
I2004t IPrebel et all 12007 al) . Additionally, when considering 
the [Z/H] MDF, we need to be mindful of issues including 
(a) giants, in general, offer a larger suite of measurable el- 
ements than dwarfs, (b) for a fixed abundance, the lines in 
giants are generally stronger than in dwarfs, thereby enabling 
measurements in giants, rather than limits for dwarfs, in many 
cases, and (c) the highest values of Z in Figure|5]; likely suffer 
from large incompleteness. Furthermore, we note that inclu- 
sion of C, N, and O abundances may considerably alter the 
[Z/H] MDF compared to our present distribution. (We empha- 
size again that throughout the present paper the MDF refers to 
the [Fe/H] distribution function unless specified otherwise.) 
Armed with sufficient observational data, MDFs can of 



course be constructed for a range of specific elements, e.g., 
[O/H] and [C/H], rather than just for [Fe/H] or [Z/H]. Such 
element-specific MDFs can then be compared with the out- 
puts of various chemical evolution models, as we did for 
[Fe/H] in Figure [3] Doing so may provide valuable insights 
into the triumphs and deficiencies of those models, and indi- 
cate ways in which they can be improved. 

3.3. The Fraction of Carbon-Enhanced Metal-Poor { CEMP) 

Stars 

In Figure [6j we again plot the raw MDF (using generalized 
histograms), but on this occasion we also include in the figure 
the MDF when restricted to C EMP objec t s, whe re we have 
used the CEMP definition of lAokiet ail d2007h ([C/Fe] > 
+0.70, for log(L/L ) < 2.3 and [C/Fe] > +3.0 - log(L/L w ) 
for log(L/L Q ) > 2.3; as opposed to the iBeers & Christliebl 
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produced directly from this panel. In panels (f) through (i), we overplot our corrected MDF adopting the error bars from the raw MDF (we do not show error bars 
below [Fe/H] = —4.25 since they extend beyond the range of these plots). 

and metallicitjH- We also note that for stars with large [C/Fe] 
ratios (and for metal-poor stars in general), a more rigorous 
chemical abundance analysis would require, amongst other 
things, model atmospheres with appropriate CNO ab undances 
and c onsideration of 3D and/or non-LTE effects dAsplundl 
120051) . Bearing in mind these shortcomings, as well as issues 
regarding selection biases and completeness of our sample al- 
ready discussed, we now comment on the fraction of CEMP 
stars. 

19 Had we updated the [C/Fe] ratio via [C/Fe]N cw = [C/Fe] Literature — 
([Fe/H]7hi s stm i y — [Fe/H]Li KratlIre ), the numbers of CEMP objects would 



( 120051) definition of [C/Fe] > +1.0). In panel (c) we show the 
percentage of CEMP stars as a function of [Fe/H], which we 
obtain by dividing the CEMP MDF by the MDF containing 
only those stars with C-measurements or C-limits below the 
CEMP threshold. ( Here we present resul ts using the CEMP 
defini tions of both lAoki et al] 120071 and iBeers & Christliebl 
|2005[ ) Using Monte Carlo simulations, as described earlier, 
we estimate the fractional uncertainty in the CEMP MDF, and 
therefore the uncertainty in the CEMP percentage at a given 
[Fe/H]. Note that for our 38 program stars from Paper I, C 
abundances (or limits) were measured from the spectra. For 
the literature sample, we were unable to conduct the necessary 
spectrum synthesis re-analysis (since we did not have access 
to the spectra), and we chose not to make any adjustments 
to these abundances based on our adopted stellar parameters 



cha nge from 16 to 19 and from 22 to 28 for the [Beers & Christlietj <2005l) 
and lAok i et al. 1 20(53) definitions, respectively. However, we note that this 
approach only includes changes to the metallicity, and does not address any 
changes in the C abundance. 
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We find a CEMP f raction of 32 ± 8% (22 of 69) adopting 
the lAoki et al.l (120071) criteria^! and 23 ± 6% (16 of 71) using 
the lBeers & Christliebl ^OOSlT criterion for [Fe/H] < -3.0. (As 
noted above, in determining the CEMP fraction we only adopt 
stars for which we had C measurements or C limits below the 
CEMP thres hold. Thus, the to tal number of stars is not 86. 
Adopting the lAoki et~aT1 (120071) criteria, we find CEMP frac- 
tions of 25 ± 8% (11 of 44) and 29 ± 15% (5 of 17) in the 
metallicity ranges -3.5 < [Fe/H] < -3.0 and -4.0 < [Fe/H] 
< -3.5, respectively. Previous est imates of the CEMP frac- 
tion below [Fe/H] = -2.0, using the lBeers & Christliebl (120051) 
[C/Fe] > +1.0 criterio n inclu de 14 ± 4% ICohen et al.1 (120051) 
9 ± 2 % iFrebel et al.1 (12006b . and 21 ± 2% iLucatello et all 
(120061) . all of which are probably comparable with our value, 
given the differences in [Fe/H] ranges for the samples. For 
the 38 program stars of Paper I, there was a bias towards 
CEMP objects. Our somewhat subjective observing criteria 
at the Keck and Magellan telescopes, as applied to an evolv- 
ing candidate list, was to (i) observe the most metal-poor can- 
didates available, (ii) in the event of similar metallicity es- 
timates, prefer giants over dwarfs, and (iii) for more metal- 
rich candidates, observe objects with prominent G-bands in 
their medium-resolution spectra, with the expectation that a 
small fraction might be C-rich, r-process enhanced stars sim- 

20 If we had considered stars with [Fe/H] < -2.80, an arbitrarily chosen 
more metal-rich boundary, we would have o btained a CEMP fraction of 29 
± 6 % (28 of 98), using the lAoki et alJ2007l definition. 



ilar to CS 22892-052, some of which might have measurable 
Th and U for cosmo-chronometric age determinations (e.g., 
iBarklem et al.ll200H ISneden et al.ll2008l) . 

Within our sample, the CEMP fraction is higher for dwarfs 
(50 ± 31%; 4 of 8) than for giants (39 ± 11%; 18 of 46). 
This discrepancy may reflect the fact that, for a fixed metal- 
licity and [C/Fe] abundance ratio, the CH molecular lines are 
stronger, and therefore more likely to yield a measurement, in 
giants than in dwarfs. That is, some of our dwarfs have such 
high [C/Fe] limits that they may indeed have [C/Fe] > +0.7, 
and thus the CEMP fraction for dwarf stars is very likely an 
upper limit. Indeed, some 23 of 3 1 (74 ± 20%) dwarf stars 
have C limits (or no measurements), compared with only 9 of 
55 (16 ± 6%) giant stars. 

There are previous reports in the literature t hat the CEMP 
fracti on rises with decreasing metallicit y (see | Carollo et ail 
120 121 for a full description). Including the lCaffau et alj d201 lb 
object, three of the four stars with [Fe/H] < -4.5 are CEMP 
objects. For our sample, of the 65 stars with -4.3 < [Fe/H] < 
-3.0 and [C/Fe] measurements, 19 are CEMP objects. Adopt- 
ing this CEMP fraction of 0.29, the probability of having three 
CEMP objects in a sample of four stars, as is the case for 
[Fe/H] < -4.5, is 0.076. While further data are clearly nec- 
essary to settle the issue, relative carbon richness at the low- 
est values of [ Fe/H] seems ubiquitous. We refer the reader to 
iCarollo et al.l (120121 and references therein), who demonstrate 
that the CEMP fraction increases from 0.05 to 0.26 ±0.03 
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as [Fe/H] as the metallicity decreases from [Fe/H] = —1.5 to 
[Fe/H] = -2.8, based on a large sample o f calibration stars 
from the Sloan D igital Sky Survey (SPSS: lYork et al.ll2000t 
iGunn et al.ll2006h . 

The above comments notwithstanding, in Figure [7] we plot 
the CEMP fraction as a function of [Fe/H] (upper panel) and 
[Z/H] (lower panel). For -4.5 < [Fe/H] < -3.0, we have three 
bins with roughly equal numbers, while the fourth bin, [Fe/H] 
< -4.5, has only 3 stars. There is no significant correlation 
between the CEMP fraction in each bin at the median [Fe/H] 
of each bin; Figure [TJa) suggests a slope of - 0.24 ± 0.22. 
Had w e included the C-normal ultra metal-poor Caffau et al. 
( 1201 11) star, we would have obtained a slope of -0.20 ± 0.19. 
For [Z/H], we use four bins with equal numbers of stars. We 
again measure the linear fit to the CEMP fraction at the me- 
dian [Z/H] of each bin. In this case, there is no significant 
correlation between the CEMP fraction in each bin at the me- 
dian [Z/H] of each bin; Figure |7Jb) suggests a slope of 0.03 
± 0.10. An important consideration is that the sample was 
selected to have low metallicity such that the stars with the 
highest Z tend to have high C abundances. Such a bias may 
potentially explain the positive trend we find between CEMP 
fraction and [Z/H]. Thus, we reiterate the need to measure O 
and N when possible to better define the metallicity, Z. Never- 
theless, we caution that the behavior of the CEMP fraction at 
lowest metallicity likely depends on the adopted "metallicity" 



definition. 



4. CONCLUDING REMARKS 



We have conducted a homogeneous abundance analysis of 
extremely metal-poor stars from an equivalent-width analy- 
sis based on high-resolution, high-S/N spectra. Our sample 
contains 86 objects with [Fe/H] < -3.0, including 32 below 
[Fe/H] = -3.5. While the completeness functions for our ~ 
30 HES program stars and the ~ 50 HK stars in the extended 
sample are well understood, the selection biases for the re- 
maining literature sample are poorly known. Nevertheless, 
our results provide an important new view of the MDF and 
CEMP fraction at lowest metallicity. 

The raw and corrected MDFs do not show evidence for 
a cutoff at [Fe/H] = -3.6. Both MDFs appear to decrease 
smoothly down to at least [Fe/H] = -4. 1 . Four stars with much 
lower metallicity are also known, three of which are present 
in our sample (th e fourth being SDSS J102915+172927; 
ICaffau et alj|20llb . 

The fraction of CEMP stars in our sample below [Fe/H] 
= -3.0 is 23 ± 6% and 32 ± 8%, w hen adopting the 
iBeers & Christliebl (120051) and lAoki et all (120071) definitions, 
respectively. The former valu e is in good agreement w ith 
previous estimates (based on the lBeers & Christliebll2005l cri- 
terion). It is unclear whether the CEMP fraction increases 
with decreasing metallicity below [Fe/H] = -3, as the appar- 
ent trend is not statistically significant (< 1 er) in the present 
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FIG. 7. — CEMP fraction versus [Fe/H] (upper) and [Z/H] (lower). In the upper panel, the lowest metallicity bin covers [Fe/H] < -4.5 while the three remaining 
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dataset. 

This study has pushed the boundary for any possible cutoff 
of the MDF down to at least [Fe/H] < -4.1, but stars below 
this metallicity are already known. Exploring the regime be- 
low [Fe/H] = -4 requires still larger samples of metal-poor 
stars, coupled with a more rigorous analysis that includes 
non-LTE effects, 3D hydrodynamical model atmospheres and, 
due to the prevalence of CEMP stars, appropriate CNO abun- 
dances. Upcoming surveys will hopefully produce significant 
numbers of metal-poor stars in the near future to address this 
need. 
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